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ABSTRACT 

The heating and cooling of transiently accreting neutron stars provides a powerful probe of the 
structure and composition of their crust. Observations of superbursts and crust cooling of accretion- 
heated neutron stars require more heat release than is accounted for in current models. Obtaining firm 
constraints on the depth and magnitude of this extra heat is challenging and therefore its origin remains 
uncertain. We report on Swift and XMM-Newton observations of the transient neutron star low-mass 
X-ray binary XTE J1709-267, which were made in 2012 September-October when it transitioned to 
quiescence after a ~10-week long accretion outburst. Within one week after accretion ended, the source 
is detected with XMM-Newton at a 0.5-10 keV luminosity of L x ^ 2 x 10 34 (D/8.5 kpc) 2 erg s" 1 . 
The X-ray spectrum consists of a thermal component that fits to a neutron star atmosphere model 
and a non-thermal emission tail, which each contribute ~50% to the total emission. The neutron 
star temperature decreases from ~160 to ~149 eV during the ~8-hour long observation, while the 
non-thermal component remains constant. We interpret this as cooling of a layer located at a column 
density of y ~ 5 x 10 12 g cm~ 2 (~50 m inside the neutron star), which is just below the ignition depth 
of superbursts. We constrain the heat generation in the layers on top to be ~0. 09-0. 16 MeV per 
accreted nucleon. The magnitude and depth rule out electron captures and nuclear fusion reactions 
as the heat source, but it can be accounted for by chemical separation of light and heavy nuclei. 
Subject headings: accretion, accretion disks - dense matter - stars: neutron - X-rays: binaries - 
X-rays: individual (XTE J1709-267) 
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1. INTRODUCTION 

Many observable properties of neutron stars are set by 
the structure and composition of their crust. Low- mass 
X-ray binaries (LMXBs) offer powerful means to investi- 
gate the properties of these layers. In LMXBs, a neutron 
star accretes matter from a late-type companion star that 
overflows its Roche lobe. This typically generates a 2- 
10 keV luminosity of Lx — 10 35-38 erg s _1 . Transient 
LMXBs only occasionally experience outbursts of accre- 
tion, while they spend the majority of their time in a 
quiescent phase with little or no matter being accreted 
onto the neutron star. The quiescent X-ray emission is 
therefore much dimmer; Lx — 10 31-34 erg s^ 1 . 

Matter that falls onto the surface of a neutron star is 
compressed to higher densities as it becomes buried by 
freshly accreted material. As a result, the matter under- 
goes compositional changes that are driven by nuclear 
reactions. These involve thermonuclear burning near the 
stellar surface, electron captures and neutron emissions 
that occur a few tens of meters inside the neutron star, 
and nuclear fusion reac tions that take place at several 
hundred meters depth ( Sato Il979t iWallace fe Wooslevl 



[lffD iHaensel fe Zduniklfl990l: ISchatz et al.H200lD . 

Initial calculations suggested that this chain of nuclear 
processes deposits an energy of ~1-1. 5 MeV per accreted 
nucleon in the neutron star crust (jHaensel fe Zdunikl 
119901 [20031 ) . However, observations of superbursts (ther- 
monuclear explosions caused by unstable carbon burn- 
ing) and the thermal evolution of neutron stars after 
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an accretion phase, both require a h otter crust than 
these heating rates can account for dCumming et aLl 
120061 : IBrown fe Cummingl l2009t IDegenaar et all 1201 lal ). 
The discrepancy may be resolved if there is an addi- 
tional source of heat from nuclear reactions or other 
processes in the outer crustal layers flGupta et al.l [20071 : 
iHorowitz et al.ll2008l: iMedin fe Cumminell201lD . " 

During outburst, the energy that is generated in the 
crust creates a temperature profile that depends on the 
depth and magnitude of the heat sources. During qui- 
escence, the gained heat is thermally conducted towards 
the surface and the core, until the crust has cooled and 
thermal equilibrium is restored. The cooling of neutron 
star crusts can be observed as a decrease in quiescent 
therm al X-ray emission after the end of an accretion out- 
burst (IWiinands et all I2IM 12001 JCackett et all [2fM 
20101: IDegenaar et all 120091 l2011a llbl: iDfaz T rigo et all 
20TTt iFridriksson et al.ll2010L I2011D . 



Comparing crust cooling observations with neutron 
star thermal evolution models yields im portant insight 
into t he properties of the cru st and core (R utlcdge et all 
2001 iShternin et all [20071 IBrown fe Cummind 120091: 



Degenaar et all 1201 latlPage fe Reddvll2012D . However, 



the current data provide only weak constraints on any 
additional heat sources. Progress can be made by observ- 
ing the temperature evolution shortly (within ~2 weeks) 
after the end of an outburst, since that can provide 
a direct measure of the dep th and magnitude of the 
heat release in the oute r crust ()Brown fe Cummingi r2009: 
IDegenaar et al.ll2011al) . 

XTE J1709-267 (1RXS J170930.2-263927) is a neu- 
tron star LMXB that was first reported in outburst 
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TABLE 1 
X-Ray Observation Log 



Observatory 


Date 


Obs ID 


Exposure time 
(ks) 


Swift 


2012 Sep 18 


35717011 


0.9 


Swift 


2012 Sep 23 


35717012 


0.8 


Swift 


2012 Sep 28 


35717013 


0.9 


XMM-Newton 


2012 Sep 30 


0700381401 


31 


Swift 


2012 Oct 16 


35717014 


0.9 


Swift 


2012 Oct 20 


35717015 


1.1 


Swift 


2012 Oct 23 


35717016 


0.9 



Notes. After excluding a short episode of background flaring, 
the net exposure time of the XMM-Newton observation was 
~30 ks for the PN and ~28 ks for the MOS1. 



in 1997 (jMarshall et all 119971 ). The source is thought 
to be associated wit h the globular cluster NGC 6293 
(jJonker et alJ l2004bl). which wo u ld place it at adis- 
tance of -8.5 kpc (jrlarrisl 11995 ILee fc Carnevl I2006T) . 
In the past 15 years, it has exhibited several accre- 
tion outbursts with a typical 2-10 ke V intensity of 
L x 2 x 10 37 (P/8.5 k pc) 2 ergs" 1 (iMarshall et al " 
1997; ICocchi et all fl99l 



Jonker et all I2003L I2004all 



Markwardt k Swankl2'004l : lRemillard2007l:iNegoro et all 
2010t lYamauchill2010l) . ~ 

The most recent outburst occurred in 2012 June- 
September (|Sanchez-Fernandez et all 120121) . Figure Q] 
displays the light curve as observed with the MA XI 
all-sky X-ray monitor. The source was active for 
~10 weeks with an average 2-10 keV flux of Fx — 
2 x 10 -9 erg cm~ 2 s _1 . Its intensity remained below the 
detection limit of MAXI from September 11 onwards, 
marking its return to quiescence (Figure [1} . 

In this work we report on a series of Swift Target of 
Opportunity observations and an XMM-Newton Direc- 
tor's Discretionary Time observation of XTE J1709-267. 
The data were obtained in 2012 September-October, and 
aimed to search for signs of an accretion- heated crust. 

2. OBSERVATIONS, DATA ANALYSIS AND RESULTS 

The observations discussed in this work are listed in 
Table [TJ All spectral data was fitted be tween 0.5-10 
keV using XSpec (ver. 12.7: IArnaudlll996[) . Interstellar 
absorpt ion was taken into account by using the tbabs 
model (jWilms et all I2000D . Throughout this work we 
assumed a source distance of D — 8.5 kpc. Quoted errors 
refer to 90% confidence intervals. 

2.1. Swift 

XTE J1709-267 was observed with Swift six times (Ta- 
ble [T]). Data reduction and analysis was carried out 
within heasoft (ver. 6.11). We re-processed the raw 
data obtained with the X-ray Telescope (XRT) using the 
xrtpipeline. Spectra and light curves were extracted 
using XSelect. 

The first two XRT observations were affected by pile- 
up. We therefore extracted source events using an an- 
nular region with inner and outer radii of 12" and 60" 
for observation 35717011, and 10" and 60" for 35717012. 
For the other four observations we used a circular region 
with a radius of 15". Background events were collected 
from a source- free circular region with a radius of 45". 



2.1.1. X-ray Spectra and Light Curve 

The first series of Swift observations were obtained on 
September 18, 23 and 28. We fitted the XRT spectra of 
these three data sets simultaneously to a simple absorbed 
power-law model. This resulted in a joint value of Ah — 
(3.4 ± 0.5) x 10 21 cm~ 2 and power-law indices of T = 
2.0 ± 0.1, 2.4 ± 0.2, and 2.6 ± 0.4 for the first, second 
and third observation, respectively (yielding x 2 = 0.93 
for 118 dof). The corresponding unabsorbed 0.5-10 keV 
fluxes are F x = (3.4 ± 0.2) x lO" 10 , (1.7 ±0.2) x lO" 10 
and (8.4 ± 1.3) x 10~ 12 erg cm~ 2 s^ 1 (see FigureP. 

The intensity and spectral shape suggest that the 
source was still in outburst during the first two Swift ob- 
servations, but that accretion may had ceased by the time 
of the third. We therefore attempted to fit this spectrum 
with the model inferred from the XMM-Newton data (see 
Section 12.2.11) . Leaving only the neutron star tempera- 
ture and power-law normalization as free fit parameters 
results in kT°° > 191 eV. Albeit high, this is not un- 
physical and may indicate that the neutron star surface 
was visible at this time. 

The second series of Swift pointings were performed 
on October 16, 20 and 23. XTE J1709-267 is not de- 
tected in these observations. We determined 95% upper 
limits on the count rate by applyi ng the prescript ion for 
small numbers of counts given bv lGehrelsl (Il986l) . Cor- 
responding upper limits on the 0.5-10 keV unabsorbed 
flux were estimated using pimms (ver. 4.6) by assum- 
ing a power-law spectrum with Ah = 4 x 10 21 cm~ 2 
and T = 2.9 (see Section l2T2~Tj) . This yielded F x < 
(3—7) x 10~ 13 erg cm~ 2 s -1 . Summing the three observa- 
tions reveals a small excess of photons at the source posi- 
tion for which we estimate Fx — 2 x 10~ 13 erg cm~ 2 s _1 . 
This is a factor of ~2 higher than the quiescent level ob- 
served with Chandra in 2002 (see Section r2.2.1l) . FigureQ] 
displays the flux light curve of the Swift observations. 

2.2. XMM-Newton 

XMM-Newton observed XTE J1709-267 for 31 ks 
starting on UT 2012 September 30 at 09:47. We ana- 
lyzed the data obtained with the European Photon Imag- 
ing Camera (EPIC). The PN and the MOS1 were both 
set in small window imaging mode. The source was not 
detected with the MOS2, which was operated in timing 
mode. The thin optical blocking filter was used. Data 
reduction and analysis was carried out using the Science 
Analysis Software (SAS; ver. 11.0.0). 

A short episode of background flaring was excluded 
by selecting data with high-energy count rates of 
<0.15 counts s _1 for the PN and <0.30 counts s^ 1 for 
the MOS1. This resulted in a total net exposure time of 
~30 and ~28 ks for the PN and the MOS1, respectively. 
We extracted source events from a circular region with a 
radius of 30", and used a source- free circular region with 
a radius of 60" for the background. 

XTE J 1709-267 is detected at an average net count 
rate of 0.568 ± 0.006 counts s" 1 with the PN and at 
0.195 ± 0.003 counts s" 1 with the MOS1. Background- 
corrected light curves were created using the tasks EVS- 
ELECT and LCCORR. The X-ray spectra and associated 
response files were extracted using the task especget. 
The spectral data were grouped to contain > 20 photons 
per bin. 
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Fig. 1. — Top: One-day averaged MAXI light curve of the 2012 outburst (2-20 keV). Bottom: Evolution of the unabsorbed 0.5-10 keV 
flux after the source intensity dropped below the detection limit of MAXI (dashed line). Squares and upper limits represent Swift data, 
and the triangle indicates the XMM-Newton observation. The dotted line corresponds to the flux measured with Chandra in 2003. 

The composite model provides a good fit, yielding 
A H = (2.9 ± 0.1) x 10 21 cm" 2 , T = 1.3 ± 0.2, and 
kT°° = 155.7 ± 2.3 eV (x 2 = 1.04 for 505 dof). The 
spectral data and model fit are shown in Figure [2] The 
inferred 0.5-10 keV unabsorbed flux is Fx = (2.33 ± 
0.02) x 10~ 12 erg cm~ 2 s _1 , which translates into a lumi- 
nosity of L x = (2.01±0.01) x 10 34 (L>/8.5 kpc) 2 erg s" 1 . 
The power-law component contributes 46 ± 3% to the 
total unabsorbed flux. By setting the power-law normal- 
ization to zero and extrapolating the NSATMOS model fit 
to the 0.01-100 keV range, we estimate a thermal bolo- 
metric flux of Fboi = (1.5 ± 0.1) x 10~ 12 erg cm~ 2 s _1 . 
The spectral shape and intensity suggest that the source 
was in quiescence during the XMM-Newton observation. 

We compare these results with the quiescent proper- 
ties measured with Chandra on 2003 Ma y 12 (Obs ID 
3507; for details, see Uonker et al.1 r2004b| ) . We reduced 
the data using the CIAO tools (ver. 4.4) and extracted 
a background-corrected spectrum. We fitted this with a 
combined power law and neutron star atmosphere model, 
with A H = 2.9 x 10 21 cm" 2 and T = 1.3 fixed to the val- 
ues inferred from the XMM-Newton data. This yielded 
kT°° = 88.6 ± 8.7 eV and a 0.5-10 keV unabsorbed flux 
of F x = (1.32 ± 0.09) x 10~ 13 erg cur 2 s~\ of which 
25 ± 10% can be attributed to the power-law emission 
(x 2 = 0.84 for 14 dof). The flux measured during the 
XMM-Newton observation was thus a factor >10 higher 
than previously seen with Chandra. 
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Fig. 2.— PN (top, black) and MOS1 (bottom, red) spectra. The 
solid line represents the best fit to a combined neutron star atmo- 
sphere (dotted lines) and power law (dashed-dotted curves) model. 
The spectra were re-binned for representation purposes. 

2.2.1. Average X-ray Spectrum 

We fitted the PN and MOS1 spectral data simulta- 
neously with all parameters tied between the two in- 
struments. A simple absorbed power-law model yields 
iV H = (4.9±0.2) x 10 21 cm" 2 and T = 2.9±0.1 (x 2 = 1.32 
for 506 dof). The fit can be significantly improved by 
adding a thermal componen t. We use a neutr on star 
atmosphere model NSATMOS (|Heinke et al.ll2007l) . which 
consists of five parameters: the temperature, mass and 
radius of the neutron star, the source distance, and a nor- 
malization that reflects the fraction of the surface that 
is radiating. We fix the neutron star mass and radius to 
canonical values of M = 1.4 M© and R — 10 km, adopt 
a distance of D — 8.5 kpc, and set the normalization 
constant at 1. As such, the only free fit parameter is the 
neutron star temperature. 



2.2.2. A Decay in Thermal Emission 

The PN and MOS1 light curves reveal that the inten- 
sity of XTE J1709-267 decreased during the observation. 
The co-added light curve is shown in Figure [3l To inves- 
tigate the cause of this decay, we divided the data into 
three equal parts of ~ 10 ks and analyzed the individual 
spectra. We used the combined power law and neutron 
star atmosphere model described in Section 12.2.11 and 
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Fig. 3. — Top: Combined PN/MOS1 light curve using a bin time of 1000 s. Bottom: Evolution in neutron star temperature along with 
a fit to a power-law decay with an index of -0.04 (dashed curve) and an exponential decay of 3.5 days (solid line). 



required the hydrogen column density to be the same 
between the three data segments (for a justification, sec 
iMiller et al~ll2009H . As illustrated in Figure IH the power- 
law index and normalization did not change during the 
observation, whereas the neutron star temperature de- 
creased. 

To better constrain the change in temperature, we re- 
peated the spectral fits with the hydrogen column den- 
sity, power-law index and normalization fixed at the val- 
ues found from fitting the average XMM-Newton spec- 
trum (Section This yielded kT°° = 1 59.8 ± 1.3, 
154.6 ± 1.5, and 149.4 ± 1.4 eV (xl = 1-05 for 450 dof). 
The temperature evolution is plotted in Figure [3J 

3. DISCUSSION 

We report on Swift and XMM-Newton observations 
obtained during the decay of the 2012 outburst of 
XTE J1709-267. The data covers a time span of 35 
days, during which the source intensity decreases from 
Lx^ 3 x 10 36 to ~ 2 x 10 33 (D/8.5 kpc) 2 erg s^ 1 (0.5- 
10 keV). The spectral shapes and measured intensities 
suggest that the first set of Swift observations trace the 
final stage of the accretion outburst, whereas the XMM- 
Newton observation and the second set of Swift data were 
obtained when accretion had ceased. 

The spectral data can be described by an absorbed 
power law that softens from T = 2.0 to 2.6 in the first set 
of Swift observations, and further to T = 2.9 when fitting 
the XMM-Newton data with the same model. Similar 
softening has been observed in other transient LMXBs 
for which high quality spectral data was obtained dur- 
ing the transition fro m outburst to quiescence (for an i n 
depth discussion; see lArmas Padilla et all [20lTl l2012f ). 
The XMM-Newton data is much better described by 
adding a thermal component, which suggests that at least 
part of the observed softening is not caused by the power- 
law emission itself, but by a soft component becoming 
apparent in the spectrum. 
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Fig. 4. — Results of time-resolved spectroscopy of the PN and 
MOS1 data. From top to bottom: the 0.5-10 keV total model 
flux, neutron star temperature, power-law index, and power-law 
normalization. 

XTE J1709-267 is detected with XMM-Newton 
at a 0.5-10 keV luminosity of Lx — 2 x 
10 34 (D/8.5 kpc) 2 erg s" 1 . The spectral data is best 
fit with a composite model consisting of a thermal emis- 
sion component that can be described by a neutron star 
atmosphere with a temperature of kT°° ~ 156 eV, and 
a non-thermal emission tail that fits to a power law with 
an index of T ~ 1.3 and contributes ~50% to the total 
unabsorbed 0.5-10 keV flux. Such a spectral shape is 
typical for neutron star LMXBs in quiescence. The Swift 
data suggests that the XMM-Newton observation took 
place within 2-7 days after the accretion ended. 



XTE J1709-267 in quiescence 
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The 0.5-10 keV flux measured during the XMM- 
Newton observation is a factor ~ 20 higher than seen 
with Chandra in 2003 (~ 1 year after the end of its 2002 
outburst), and the inferred neutron star temperature is 
higher by a factor of ~ 1.7. The subsequent Swift obser- 
vations show that the flux had decreased, but was still 
a factor of ~2 above the level detected with Chandra. 
This suggests that (part of) the crust became substan- 
tially heated during the recent ~ 10- week outburst, and 
was cooling in quiescence. Observations within the next 
year(s) may detect a fur ther decrease in thermal emission 
l|Degenaar et al.H2011al) . 



3.1. Constraining the Heat Release in the Outer Crust 

Most remarkably, we found that the intensity of 
XTE J1709-267 decreased during the ~ 8-hour long 
XMM-Newton observation. Time-resolved spectroscopy 
suggests that the thermal emission decreased, while the 
non-thermal component remained constant. This shows 
that the neutron star cooled during the observation. We 
use the neutron star ther mal evolution calculations of 
iBrown fe Cummi ng (2009) to interpret this result. 

At any given depth inside the neutron star it takes 
a certain characteristic time for heat to diffuse to the 
surface, i.e., for the layer to cool. The observed de- 
crease in neutron star temperature from kT°° ~ 160 
to ~ 149 eV can be fitted with an exponential decay 
ti me of ~ 3.5 ± 0.7 days (F igure EH ■ Using equation 9 
of lBrown fe Cummind (2009), we find that is the charac- 
teristic cooling time scale for a layer located at a column 
depth of y ~ 5 x 10 12 g cm -2 (corresponding to a matter 
density of p ~ 2 x 10 9 g cm~ 3 , or a depth of ~50 m inside 
the neutron star). This is just b elow the depth where s u- 
perbursts are thought to ignite (jCumming et al.| [2006). 

The decrease in temperature can also be described by a 
power-law decay with a slope of —0.04 ±0.01 (Figure [3]). 
This can be used to dete rmine the heat flux i n the c rust; 
using equations 8-11 of IBrown fe Cummind (|2009f ). we 
find -Fcrust = (1.5-2.5) x 10 21 erg cm -2 s -1 . Combining 
this with the mass-accretion rate provides an estimate of 
the energy that was released in the crust during outburst. 

The average intensity observed with MAXI during 
the 2012 outburst translates into a 2-10 keV luminos- 
ity of Lx ^ 2 x 10 37 (D/8.5 kpc) 2 ergs" 1 . The 
bolometric luminosity of neutron star LMXBs is typ- 
ically a factor of ~ 2.5 higher than m easured in the 
2-10 keV band (|in 't Zand et all 120071 ). which would 
imply an average accretion luminosity of L^ cc ~ 4 x 
10 37 (D/8.5 kpc) 2 erg s _1 . This gives an estimated mass- 
accretion rate of M oh = RL acc /GM ~ 2.1 x 10 17 g s _1 
(for R = 10 km and M = 1.4 M Q ). Assuming that the 
emission is isotropic, the mass-accretion rate per surface 
area is rh b — 1.7 x 10 4 g cm~ 2 s _1 . We then find a 



heat release of F CTUBt /rh b — (9 — 15) x 10 16 erg g" 
0.09 — 0.16 MeV per accreted nucleon. This heat must 
have been released in the layers above the cooling layer 
that is probed by our observation, i.e., at a depth of 



i/<5x 10 12 g cm~ 2 . 

We can compare these constraints with the differ- 
ent mechanisms that have been put forward as possi- 
ble sources of extra heat release in the neutron star 
crust. The shallow depth (y < 5 x 10 12 g cm~ 2 ) rules 
out nuclear fusion reactions, since the se occur much 
deepe r in the crust (at y > 10 14 g cm" 2 ; iHorowitz et al.l 
2008). Electron captures do occur at the inferred 
depth, but these are expected to generate a heat of 
~ 0.01 MeV/nucleon (|Gupta et all 120071 : lEstrade et all 
l20ll . which is a factor of ~ 9 — 16 lower than we 
obtain. To get a heat release consistent with electron 
captures would require that XTE J1709-267 was accret- 
ing at the Eddington limit in 2012. However, as suming 
L cdd = 3.8 x 10 38 erg s" 1 (|Kuulkers et al.l 120031) would 
place the source outside the Galaxy at a distance of 
~ 42 kpc (~ 34 kpc from the Galactic center), which 
is unlikely. 

Chemical separation in the crust provides an alterna- 
tive mechanism for the heat release. Just below the sur- 
face, the crust is composed of a variety of ions that ef- 
fectively behave as a liquid. This is commonly referred 
to as the neutron star 'ocean' and is thought to be the 
ignition site of superbursts. The matter freezes out at 
the bottom of the ocean, at y ~ 10 12 g cm -2 , becoming 
part of the solid crust. Molecular dynamics simulations 
suggest that the freezing of a mixture of elements causes 
chemical separation, with light nuclei mainly being left 
in the liquid phase and heavier nuclei p referentially be- 
comi ng incorporated in the solid crust (jHorowitz et all 
120071 ). The retention of light elements in the liquid ocean 
drives a buoyancy that can cause a convective heat flux 
of ~ 0.01 — 0.2 MeV per accreted nucleon, depending on 
the composition ( Medin fc Cummindl201ll ). 

We conclude that within our current knowledge, elec- 
tron captures and nuclear fusion reactions can not be the 
source of shallow heat that we find for XTE J1709-267. 
Instead, chemical separation of light and heavy nuclei can 
account for both the depth and the magnitude of the heat 
source that we infer from our XMM-Newton observation. 
Therefore, this process likely plays an important role in 
heating the neutron star crust. It can potentially supply 
the extra heat that is required to resolve the discrepancy 
between current heating models and observations of su- 
perbursts and crust cooling of accretion-heated neutron 
stars. 
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